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Abstract. In the framework of spatially averaged inhomogeneous cosmologies in 
classical general relativity, effective Einstein equations govern the dynamics of averaged 
scalar variables in a scale-dependent way. A particular cosmology may be characterized 
by a cosmic equation of state, closing the hierarchy of effective equations. In this 
context a natural candidate is provided by the Chaplygin gas, standing for a unified 
description of dark energy and dark matter. In this paper, we suppose that the 
inhomogeneous properties of matter and geometry obey the Chaplygin equation 
of state. The most extreme interpretation assumes that both dark energy and 
dark matter are not included as additional sources, but are both manifestations of 
spatial geometrical properties. This feature is an important conceptual difference in 
comparison with the standard approach of a Friedmann-Lemaitre-Robertson- Walker 
universe filled with dust and another fundamental field characterized by the Chaplygin 
equation of state. We finally discuss the consequences of the resulting scenario for 
effective cosmological parameters in order to establish the framework of a future 
confrontation with observations, and we note that the standard Chaplygin gas may 
not be ruled out by them. 



PACS numbers: 04.20.-q, 04.40.-b, 95.35. +d, 95.36,+x, 98.80.-Es, 98.80.- Jk 



1. Introduction 



Does an inhomogeneous universe evolve on average like a homogeneous solution in the 
framework of general relativity? This question is not new [26J and naturally emerges 
in view of the nonlinearity of the theory and, in particular, from the generally non- 
commuting operations of averaging and time evolution [27]. The main difficulty to 
answer it resides in the notion of averaging and in its construction (see, e.g., [19], section 
2.2 of |4j, and references therein). 

Our universe is supposed to verify the strong cosmological principle which demands 
homogeneity and isotropy at all scales. This standard approach, known as Friedmann- 
Lemaitre-Robertson-Walker (FLRW) cosmology, is widely used in order to describe 
the dynamics of our universe and the formation of its constituents. It however leaves 
in suspense an explanation about the origin of dark energy and dark matter, which 



Chaplygin Gas and Inhomogeneous universe 



2 



respectively represent in this model about 3/4 and 1/4 of the total content of the 
universe. This last point might actually reveal a symptom of a deeper problem linked 
to this approach. Indeed, in FLRW cosmology one determines background quantities 
regardless of the scale and makes them evolve according to a homogeneous-isotropic 
solution of Einstein equations. Our first query could be reformulated in order to note the 
central aspect of this issue: Are the background quantities well defined within standard 
cosmology, i.e. as a suitable average over the inhomogeneities? Is their evolution well 
approximated in this framework, i.e. is the time dependence of the homogeneous- 
isotropic averaged state well approximated by a homogeneous-isotropic solution? 

We shall adopt an approach that averages, in a domain-dependent way, the scalar 
parts of Einstein's equations with respect to synchronous free-falling observers in a 
dust model [T3l[T4]. a realization of the averaging problem that does not answer the 
above questions in the affirmative. The average evolution of an inhomogeneous universe 
differs from the evolution of a homogeneous one; in other words, even if we are entitled 
to describe structure formation in terms of perturbations of a background, this latter is 
generally not a member of the homogeneous solutions (see also [35l|36]). This difference 
of evolution is driven by the non-trivial geometrical structure of an inhomogeneous 
space, featuring deviations that are known as "backreaction" . These backreaction effects 
can act on average, at least qualitatively, as the dark components. 

The set of equations obtained within this approach should be closed to derive the 
evolution of all the involved quantities, namely the effective scale factor, the averaged 
scalar curvature deviation and the kinematical backreaction variable. In recent papers, 
attention was turned to a closure under the assumption of global constraints such as 
a globally stationary universe OdE], or by exploring the solution space with exact 
scaling laws for the backreaction and the averaged scalar curvature [HI EE], or by 
symmetry requirements such as spherical symmetry (e.g. [IT], [21], [281 E3 EH HH] , and 
references therein). In this work we want to choose the closure relation by focusing 
on the particularity of the model to unify the dark components through backreaction. 
According to this point of view, the Chaplygin gas (CG) seems to be an interesting lead 
since it unifies dark matter and dark energy in only one fluid, obeying an exotic equation 
of state [3T1I32P 34] . This unification is made through the evolution of this particular fluid 
and it can be extended to a unification where both dark components are simultaneously 
modelled thanks to the scale dependence of our approach^]- These points motivate us 
to build a model in which the generically existing coupling between the backreaction 
and the averaged scalar curvature deviation, which encodes the particular geometrical 
structure evolution of an inhomogeneous universe, is furnished by a scale-dependent CG 
equation of statjjjl 

| A 'simultaneous' unification has also been proposed in the context of an inhomogeneous, but 
fundamental CG in [5]. 

§ A study of inhomogeneous spherically symmetric spacetimes, presenting nonlinear perturbations 
constructed from the fluctuations of local variables with respect to background quantities called quasi- 
local scalars, has been given in [33]. As an example the CG is employed in [33] and relates the quasi-local 
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In section [2] we introduce the equations that govern the average evolution of an 
inhomogenous universe model and briefly discuss them. We present in section [3] the 
CG and some of its properties. In section H] we suppose that the backreaction and the 
averaged scalar curvature deviation are coupled according to the CG equation. The CG 
then describes a particular geometrical structure of the inhomogeneous universe and 
does not correspond to any fundamental field. Finally, in section |5l we reformulate the 
results obtained in terms of effective cosmological parameters; we compare this model 
to the Friedmannian framework and we study in particular the acceleration of a spatial 
domain. 



2. Effective description of inhomogeneous universe models 

Restricting attention to a universe filled with irrotational dust, i.e. irrotational 
pressureless matter, we spatially average the scalar parts of Einstein equations (the Ha- 
miltonian constraint, Raychaudhuri's equation and the continuity equation) with respect 
to a collection of comoving (generalized fundamental) observers over a compact, restmass 
preserving spatial domain V, and obtain the following set of equations ( [T3l|lll[T9] > IS] 
for details): 

a v \ 2 8itG _ (R) v + Qy m 



a v J 3 x lu 6 

^ + ^(,) c = %, (2) 

(qYv + Z — {q) v = °> (3) 
-4- (QvalY + \ ({R) v a 2 v y = 0, (4) 
where a v is the effective volume scale factor 



v. 



with Vd 1 the initial volume of the domain and Vv{t) its volume at a proper time t, 
(g) v = M a^ 3 /V-Di the density of irrotational dust averaged over V, (R) v the spatial 
scalar curvature averaged over T> and Qx> the kinematical backreaction 

Qv{t) := \ {{9 - (6) v )% - 2 (a% , (6) 



with 9 the rate of expansion and a := y \° l ^ i3 the rate of shear with the shear tensor 
components a^. 

Equations p]) and ^ govern the kinematics of the effective scale factor and 
equations ([3]) and (jlj) express the conservation law for the dust matter and the 

variables' pressure and density. The reader may find connections between quasi-local variables and our 
averaging procedure in }43| . 
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backreaction terms, respectively. It is important to point out that (R) v might evolve 
differently from its Friedmannian counterpart (as we shall discuss in subsection 14. 3p . 
Upon introducing the averaged scalar curvature deviation Wd := (R) v — Gk^a^ 2 , we 
may rewrite equation ([1]) in the form 

av\ 2 , StiG Wt> + Qv 

- v ) + ^ "^ ( ^ = — • (?) 

Looking now at equations ([2]) and (J7J one should note that both the backreaction and 
the averaged curvature, through Wp, induce a change in the averaged dynamics of the 
domain in comparison with the Friedmannian framework. Equation (j2J) states that a 
positive backreaction contributes to accelerate the expansion of the domain and then 
plays against gravity: Qd > effectively mimics a dark energy behaviour over T>. The 
domain will actually undergo an acceleration of its expansion only if the 'intensity' of 
dark energy is sufficient, which is the case when Qx> > 47rG(£>)x>. A negative backreaction 
contributes to decelerate the domain expansion and therefore adds to gravity: Qv < 
effectively mimics a dark matter behaviour over T>. For the averaged model we may 
suppose that backreaction acts as dark matter on small scales (e.g. galaxy cluster 
and void scales) and as dark energy on the largest scales (CMB and high-redshift 
supernovae). In the present work, this differentiation with respect to the spatial scale 
will however not be made explicit. An explicit multi-scale dynamics can be formulated 
to refine such a description [26]. 

We shall assimilate here the properties of the spatial geometrical structure to a 
domain-dependent CG. To this aim we first describe the backreaction variables in terms 
of an effective perfect fluid whose energy density and pressure reacfjj] 

We stress here that, since it is an effective description, this fluid does not have to satisfy 
any energy conditions (as discussed in [H]). We therefore reformulate equations ([TJ and 
(EJ) casting them into Friedmannian form 



av\ , k Vi 8%G 



. + - (te)i> + fi?) = °' (9) 
~-+^((Q)„ + e? + 3p?) =o. (10) 

Using the last two equations together with equation Q one obtains the conservation 
law for the backreaction fluid 

fi? + 3^ ( Q ?+p?)=0, (11) 

|| In this work we prefer to consider the deviation term to describe the fluid, instead of the 
full averaged scalar curvature (R)-p. First, (Qx>jWx>) incorporate the deviation from a general 
Friedmannian model, being equivalent to the pair (Qv, (R)x>) only in a zero-curvature Friedmannian 
model. Second, the kinematical backreaction and the curvature deviation both vanish on the 
background and are gauge-invariants (as shown to second-order in perturbation theory [39 ). 
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which, if written out, reflects the generic coupling between the curvature deviation and 
the kinematical backreaction — a simple reformulation of equation (j3J): 

\ (Qva%\ + \ (W v a 2 v y=0. (12) 
3. The Chaplygin gas 

The CG is a perfect fluid obeying the state equation 

A . . 

Pch = , (13) 

Qch 

where g c h > and p c h are respectively the energy density and the pressure of the fluid in 
a comoving frame and A is a positive constant. It was first introduced as a cosmological 
fluid unifying dark matter and dark energy by Kamenshchik et al [34J and has been since 
widely studied in this context (see, e.g., |2ll^UniDinil2nil221l2Hll2SU2S-[H2] ) - The equation of 
state (I13p has also raised interest in particle physics thanks to its connection with string 
theory [12] and its supersymmetric extension [33J. The generalization of the CG [5||3^]. 

Pch = , (14) 

Qch 

with a a free positive parameter, is commonly used in cosmological models; however in 
the present work we shall consider for simplicity the case a — 1, i.e. the standard CG0. 
Assuming that the gas verifies the energy conservation law over a spatial domain £>, 



— [Qch 
a?) 

we obtain, making use of relation (JT3"]) . the expressions 



eSl + 3— (el+P? h ) = o, (is) 



Qch = \ A v + — , P ch = . ===== , 16 

where Bd = gf h . — Ax> determines the initial conditions of the CG and both Ax> 
and Bx> depend on the domain. Note that equations f|T6l) describe the evolution of 
a homogeneous CG which is the one of interest in our work since the backreaction terms 
are, due to the averaging procedure, homogeneous over a spatial domain. However, 
since we average over inhomogeneities, there certainly exist interesting links to the 
inhomogeneous CG (for an investigation of the latter the reader is referred to [H]). 

We now recall briefly of the different aspects of the CG according to the sign of Bx>- 

A generalization of our ideas, using equation (|14|) . is straighforward. We emphasize that, even if the 
standard CG does not seem to well fit with observations in a FLRW model, this must not be the case in 
our approach, since observational data have to be reinterpreted before the need for such a generalization 
is justified (cf section [S]) . 
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3.1. Choosing a positive integration constant 

We follow [34] to present the behaviour of the CG for a positive B-d- Expressions (jTBl) 
become for small values of the scale factor, a v <C Bd/At>, 

r B- 



Qch~—5-> Pch~°> ( 17 ) 
a v 

which indicates that the CG can behave as a dark matter component. For large values 
of the scale factor, a v 3> Bt>/Ad, it follows that 

*&~V^, p£~-\/^, (18) 

which reflects the dark energy-like behaviour of the CG in its last stage. Finally, one 
may also develop (fl6|) for large values of a© to obtain 



e * ~ ^ + — ^a^, ^~-v / A I? + - 7 ^a^. (19) 

Between the phases (IT7|) and (ITS1) the CG can be seen as a mixture of a cosmological 
constant and a stiff fluid whose pressure and energy density are equal. 

To resume, for a positive B-d, the CG acts first as dark matter and then as dark 
energy whose state equation evolves towards the one of a cosmological constant. 

3.2. Choosing a negative integration constant 



As already noticed in [42J the CG presents another interesting feature for a negative 
Bx> since its density increases with the scale factor. In this situation it plays the role of 
phantom dark energy, and it evolves at late times towards a cosmological constant. For 
the pressure and energy density to be well defined one needs 



A v + ^ > o a% > . (20) 



It therefore exists a minimal value for the scale factor, a min = (— Bd/At,) 1 / 6 , implying 
that this case describes a bouncing universe model at early times. 



4. Evolution of the kinematical backreaction and the curvature deviation 

4-1. Backreaction fluid as a Chaplygin gas 

As we have outlined in section 121 the backreaction fluid inherits a simultaneous 
unification of the dark components thanks to its scale dependence. We want to build a 
model in which this fluid is assimilated to the Chaplygin gas in order to physically shape 
the behaviour of the backreaction and the curvature deviation on a given spatial scale, 
while also allowing for a metamorphosis of the 'dark character' through evolution. In this 
model, the Chaplygin gas emerges from the inhomogeneous structure of the universe, 
and is not related to any fundamental field. We then consider that the backreaction fluid 
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responds to the scale-dependent state equation ( TT31 over a spatial domain V. Using 
definitions (JS}, we rewrite expression f lT3|) in the form 

{Wv + Qv)(^y- - Qv) = (l6nG) 2 A v . (21) 

Since no energy condition has to be verified by the backreaction fluid, we may also 
consider the Chaplygin equation of state (fl3]) with a negative energy density (and hence 
a negative Av if one wants to preserve the negativity of the pressure). However, we 
shall restrict ourselves in this paper to the case where the backreaction fluid satisfies 
f fT3|) under the usual conditions, i.e. with a positive energy density and a positive A jj. 
We thus have to respect, in view of expression fl2T|) . the following constraints, which we 
call the Chaplygin fluid constraints: 

Wv 

Wv + Qv < ; - Q v < , (22) 

or, equivalently, 

Wv 

Wv < , — ^ < Qv < -Wv ■ (23) 

In this situation, the curvature deviation of any domain is negative at any time whatever 
its dynamics could be. We also note that Qv and Wv evolve in such a way that relation 
(|2T|) is always satisfied (see figured]). Rewriting equation ( 12T|) as a function Wv(Qv, Av) 
under the Chaplygin fluid constraints, we derive, for a given Av, the maximal value of 
the deviation term Wm = —2\-kG^A~v obtained for Qw M = SnGy/Av- We shall see 
in subsections 14.31 and 14.41 that these values form an attractor for the dynamics of the 
system. 

4-2. Exact evolution equations 

Inserting expressions (|8]) into the evolution laws ( Fl6l) results in 

^JA V + B v /a% 2<4, 
„ 8tvG , . Bv \ 

^ Av + Bv/a v ^a v 

We may express the constants Av and Bv in terms of the initial values of the 
backreaction Qi and the curvature deviation Wj over V, 

A - = ^^ + ^ {Y Y- Q ^ < 26) 

Bv = J^ffl (Wi + Qi) + ft) , (27) 

and determine the evolution of Wv and Qv in the form 

+ We do not impose here any other energy conditions on the backreaction fluid. 
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curvature deviation 



Figure 1. Each curve draws relation (j!?Tj) for a different Ad (its value arithmetically 
increases from the curves on the right to the ones on the left) under the Chaplygin 
fluid constraints. For a given Ad, whatever their evolution could be, Qd and W-p 
run on the corresponding iso-Ad curve. The dots represent the maximal value of the 
curvature deviation for each Ad- A negative Qd mimics dark matter and a positive 
Qd dark energy. 



a (3 + cry a 



v 



v 



Q 



V 



2 (a/3 + 2 cry a^ 6 ) 1/2 
1 cr/3 — cry a^ 6 



2 (a/3 + 2 cry a P 6 ) 1/2 
where the new terms are defined as 

a := W t + Qi , 



:= -^r - <?« 



(28) 
(29) 

(30) 



The Chaplygin fluid constraints imply a < and < 0. The evolution of and Wx> is 
entirely determined by the initial values Qi and Wj of the domain or, equivalently, by a 
and (3. The opposite sign of 7 gives the sign of B V) equation (|2T|) . and therefore defines 
the behaviour of the Chaplygin backreaction fluid. 



4-3. Evolution of the curvature deviation 

The rate of change dWx>/dax> shows that the curvature deviation grows with a©. For 
large values of the scale factor, a% 3> |7//3|, we get 



•|M) 1/2 



(31) 
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(32) 



where Wm has been introduced before. Thus, for a given Ap, the curvature deviation 
increases towards its attractor Wm, for which the departure from the Friedmannian 
curvature is maximum. For 7 = the curvature deviation is initially set to Wm and 
then does not evolve. It is important to remark that there exists an infinite number of 
couples (Qi, Wj) or, equivalently, (a, (3) which yield the same Ax> under the Chaplygin 
fluid constraints. For all of them W© will live on the same iso-Ap curve and will tend 
towards the same attractor (see figure [1]) . 

4-4- Evolution of the backreaction 

For large values of the scale factor, a% ^> |7//3|, we have 

^~^M) 1/2 , (33) 

or, equivalently, using equation (1261) . 

Q v ~ 8nGy^ = Q W M > , (34) 

where Qw M has been introduced before. The backreaction tends at late times towards a 
domain-dependent cosmological constant whose value is given by the initial conditions 
a and /3 on the domain. The same remark as in subsection 14.31 can be made for Qx>- 
One needs to distinguish the following cases (see figured]). 

(i) For 7 > the backreaction is always positive and hence only acts as dark energy 
over the domain. In this situation, since dQx> / da-p is negative, Qd behaves as dark 
energy whose intensity decreases until reaching the attractor Qw M . We recall in 
this case the existence of a minimal scale factor a£} in = (— 2 / y//3) 1 ^. a < implies 
a™ in < 1. 

(ii) For 7 = the system is initially set on the attractor; thus it does not evolve. The 
backreaction always acts cosmological constant. 

(iii) for 7 < we have a positive dQv/da-p- Two subcases arise according to the initial 
value Qi. 

(a) for < Qi < — Wj/3 the backreaction acts as dark energy whose intensity 
increases until reaching the attractor Qw M ; 

(b) for Wj/3 < Qi < the backreaction changes its sign during its evolution. 
It first behaves as dark matter whose intensity decreases, and then as dark 
energy whose intensity increases until reaching the attractor. The moment of 
the transition dark matter-dark energy depends on the initial values of the 
domain since a% = (7//3) 1 / 6 . 
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We conclude by emphasizing some of the characteristics of the model. First, the CG 
relates the backreaction terms of our inhomogeneous universe model and therefore 
furnishes a particular manifestation of its inhomogeneous structure. Our model does 
not suppose the existence of any other fundamental field, contrary to the approach 
of a FLRW universe filled with dust and CG. Second, the role and the evolution of 
the backreaction depend on the domain considered and on its initial values. As a 
consequence this model could be seen as an effective multi-scale model. For instance 
Qx> might act as dark matter on small scales (situation (iii)-(b)) with a dark energy 
transition occuring at very late times, and as different kinds of dark energy on larger 
scales (situations (i), (ii) and (iii)-(a)). 

4-5. Another approach: the backreaction fluid as a scalar field 

The backreaction fluid might also be described by a minimally coupled real scalar field 
0x>, called the morphon field [18], evolving in an effective potential Uv{4>v), as follows: 



+ P?:=ffi-U v , (35) 

where e = +1 for a standard scalar field (with a positive kinetic energy) and e = — 1 for 
a phantom scalar field (with a negative kinetic energy). The last expressions together 
with relations flH]) give 

The system evolves towards the maximal value of the curvature deviation Wm, as 
seen in subsection 14.31 which corresponds to the minimal value of the potential. With 
correspondence fl36l) the integrability condition ffl2]) implies that 0x>, for <fix> 7^ 0, obeys 
the scale-dependent Klein-Gordon equation 

4> v + 3— <k> + Ctt-Uv = . (37) 
a v 0<p v 

The above correspondence allows us to interpret the kinematical backreaction effects in 
terms of the properties of scalar field cosmologies, notably quintessence or phantom- 
quintessence scenarii that are here routed back to models of inhomogeneities. The 
morphon field may also be characterized by the domain-dependent equation of state 
p® = w^gV , which assumes in our model the form 

< = — z 7 -— ■ ( 38 ) 

1 + 2ja v e 

For 7 < , the morphon field acts as dark matter and 'standard' dark energy as 
mentioned in subsection 13.11 In this situation we have —1 < < and — > — 1 + at 
late times. For 7 > , it behaves as phantom dark energy as mentioned in subsection 
13.21 In this case we have, since ax> > a© 111 , < —1 and — > — 1~ at late times. We 
stress again that the phantom character is an effective property in our description; no 



Chaplygin Gas and Inhomogeneous universe 



11 



fundamental phantom field is assumed to exist [TB]. Finally, for 7 = 0, the morphon 
field mimics a scale-dependent cosmological constant and we have = — 1. 
Using equations ( |28i) and ( 1291) . relations ( 1351) become 

«4 2 - — alai ° (39) 

167rG (a/3 + 2a 7 a^ 6 ) 7 

The scalar field dynamics can be reconstructed by evaluating the following integral, e.g. 
to find Uj)((pv)- 

<Kf = -T^ ~r U 7? 7T72 , ( 41 ) 



4vrG 



Mo , -fi\i/2 -2! -fiN 1//2 

— a^, + (a/3 + 2 «7 a P j —fSkvflv + 2 cvy a v 



where the prime denotes the derivation w.r.t. the volume scale factor. This relation 
does not seem analytically integrable in the general case. In the vacuum with a zero- 
Friedmaniann curvature (k^ = 0), expression (j4"Tj) becomes 

'(^) 2 - * IT* -« • < 42 > 

47rG a/3 + 2 a7 a v 

From this last relation, we get the following expression for a standard scalar field (e = +1 
and 7 < 0): 



1 / A -+- Ba~^ 
(pv(av) = T ^24^Q arccosh y J- 2 - + ^0 , (43) 

where </>o is an integration constant. The potential is then written [M] 

M<M = ^ ( cosh v / 24^G(0 P - O ) + - 1 —) ■ (44) 

327rG \ cosh V247rG(0x, - 0o) / 

For a phantom scalar field (e = — 1 and 7 > 0) we get from the relation 



/ \ 1 A + Baj? , , 

<Ma©) = T^== arccos y + O , (45) 



and the potential reads 

Uv((h) = ^rT^V 7 ^ (cos V247rG(0p - <j> ) + 

OZ7TG \ 



(46) 



cos v / 247rG(0x> — 0C 

Note that this last expression can be obtained directly from relation ( |44p by the simple 
reparametrization </>£> — <Po ~ > i(<h> ~~ 0o)- This is obvious in view of the form of the 
kinetic term of the morphon field: for a standard scalar field it is written 0|, and for a 
phantom field — 0|>. 

In the following section, we reformulate the different results obtained in this section 
in terms of effective cosmological parameters and we study the dynamics of the model. 
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5.1. Constraints and evolution equations for the cosmological parameters 
Expressed through the domain-dependent cosmological parameters 

S%:=%£(e)v, ^=-7^V, (47) 

w 6#|,' Q ml' y ' 

where Hx> = dx>/ax> is the volume Hubble functional, the averaged Hamiltonian 
constraint (HI) assumes the form [19] 



n£ + n? + n5, + ng = i. (49) 

The Chaplygin fluid constraints ( 122]) may be expressed as 

Oft + ngX), ^-fig>0. (50) 
Equations ( 1281) and ( 129]) furnish the evolution laws for and fin : 



2 (fi a ^ + 2 Q a Q^ a v 6 ) 1/2 

where the new terms are defined by 

a (3 7 1 



The sign of f2 7 determines the behaviour of Qyy and VIq. One may reformulate 
constraints (|50|) in terms of the initial conditions, 

< f2 Q , < Sip. (54) 
5.1?. Evolution of the Hubble functional 



Equation ( 1491) together with equations ( 1511 and (1521) provide the evolution equation for 
the Hubble functional (see appendix A for the study of Hx> in a particular case): 



tt2 

-^g = Ql a^ 2 + W m + (fi Q fig + 2 fi Q fi 7 a^ 6 ) 1/2 . (55) 



if^ tends, at late times, towards 



ft ~ (0^) 1/2 • (56) 
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5. 3. Evolution of the cosmological density parameters 
One finds, for large values of the scale factor, 

Vq~-\- ( 57 ) 

It is interesting to note that the cosmological parameters tend towards a value 
independent of the initial conditions (see appendix A for the study of Qy^ and Qq in a 
particular case). We may also introduce the cosmological parameter Q x , the so-called 
X-matter, defined as 

H 2 1/2 

We have noticed in section [2] that the kinematical backreaction might explain the origin 
of dark matter and dark energy, and we have presented its behaviour in subsection 
14.41 If one wants to compare our model to a scale-dependent Friedmannian cosmology 
in terms of cosmological parameters, Q x has to be considered as the origin of the dark 
components, instead of VLq alone. This simply means that the curvature deviation, since 
it also participates in the departure from the Friedmannian framework, may also 
act, qualitatively, as the dark components. In fact we know from previous work that 
the curvature deviation is actually quantitatively more important than the backreaction 
term itself [IH] • 

We depict in figure [2] two situations to illustrate our model. The initial moment is 
chosen to be the epoch of the CMB, and we set = C0, fi^ = 1 - 1 ®~ 5 and = 10 ~ 5 - 
In the upper figure, the X-matter only stands for dark energy; Q x corresponds to Qp E 
(the dark energy in the concordance model) and to Qb + ^dm (respectively the 
baryons and the dark matter in the concordance model). In the lower figure, the X- 
matter stands for dark energy and dark matter (since it can play both roles, it is regarded 
in this case as a result of different contributions on different scales); Q x corresponds to 
+ and to fib- The two situations describe a universe initially close to a 
homogeneous-isotropic Friedmannian universe (fl x ~ 0) with zero-curvature (flj? = 0). 
This low deviation, however, becomes larger with the growth of the scale factor to reach 
in the first case f^(ai) ~ 1000) = Qx = ~ 0-^2 (where the superscript denotes the 
'today'-value), and in the second case VL x (a v ~ 1000) = Vt Q x = Vt° DE +Vt Q DM ~ 0.950. The 
origin and magnitude of the dark energy (first situation) and of both dark energy and 

* This choice allows us to roughly compare our model to the concordance model. In this situation we 
have Wv = (R) v - 

jj If we choose the initial moment to be the CMB epoch, a-p ~ 1000 is chosen to roughly correspond to 
the 'today'-value of the scale factor in a Friedmannian cosmology. Under the assumption that the metric 
of our universe does not significantly differ from a Euclidean metric, the scale factor, calculated through 
this metric (equation (0)), evolves in our model approximately like its Friedmannian counterpart. Thus, 
we are entitled to assume a^, ~ 1000. However, backreaction terms, encoding the inhomogeneities, 
involve first and second derivatives of the metric and can therefore not be neglected, even if the metrical 
amplitudes are considered to be near-Euclidean (see the discussion in |17j and numerical applications 
therein). 
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Figure 2. Evolution of the cosmological parameters Sl„ l7 f2^, fly^ and fig w. r. t. 
the scale factor. We set ttf = 0, Vt i m = 1 - 10~ 5 and fl^ = fl a = 10~ 5 where 
the initial moment is the CMB epoch. Upper figure: Qp = 5 • 10 -13 ; we have 
n%(a v ~ 1000) = n° DE ~ 0.72 (dot) and fi£(az> ~ 1000) = 0° + ft£, M - 0.28 
(square). Lower figure: Hp = 3 • 10~ n ; we have Q^f ( a i> <~ 1000) = + H° DM ~ 0.95 
(dot) and fi£(a D - 1000) = ~ 0.05 (square). 



dark matter (second situation) are then entirely explained by the particular geometrical 
structure of an inhomogeneous space obeying the Chaplygin equation of state. 
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1 10 20 30 40 50 

a D 

Figure 3. Evolution of wV w. r. t. the scale factor. Slf — 0, fl l m = 1 — 10~ 5 and 
SV X = fi Q = 10~ 5 ; the initial moment is the CMB epoch. Case (a): ftp — 5 • 10 -13 
and Vl^r — Qp E . Case (b): ftp = 3 • 10 -11 and Sl^ = ^1 1 q E + ^jjm- ^ n both situations 
quickly evolves towards —1, which corresponds to the cosmological constant -like 
behaviour of the backreaction fluid. 



5.4- Another approach: the morphon equation of state 

Using the effective cosmological parameters, we rewrite equation fl38|) as 

< = 4 • (59) 

lip 

The behaviour of is known in terms of cosmological parameters simply by replacing 
7 with — £7 7 through subsection 14.51 (see also figure [3]). 

5. 5. Dynamics of the fluid volume 

Depending on the initial conditions, the domain T> might undergo a decelerated or an 
accelerated expansion at different periods of its evolution. To learn how it behaves we 
throw an eye on the volume deceleration parameter 



to = -4^ = |^ + 2fig. (GO) 



At late times, since qv ~ — 1, the fluid volume undergoes an accelerated expansion, 
we are indeed dealing with a cosmological constant in this situation. Three different 
cases occur for the dynamics of the domain (see figure HJ and appendix A for a detailed 
analysis in a particular case): 
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case (a) 
case (b) 



1 200 400 600 800 1000 1200 1400 1600 1800 

a D 

Figure 4. Evolution of the deceleration parameter w. r. t. the scale factor, ttf = 0, 
SV m = 1 — 10~ 5 and Vt l x = Q, a = 1CP 5 ; the initial moment is the CMB epoch. Case (a): 
Sip = 5 • 1CT 13 and n% = VL^ E . Case (b): ftp = 3 • 1CT 11 and {1% = VL^ E + ft^ u . The 
evolution of the deceleration parameter is only shaped by fi^ and fig (fi^ docs not 
contribute to it, see equation (|60|) ). In both situations the expansion of the domain is 
first decelerated {qx> > 0), then accelerated (qv < 0). 



(i) its expansion is first decelerated then accelerated; 

(ii) its expansion is always accelerated; 

(iii) its expansion is accelerated, then decelerated and again accelerated. 

The latter case is an interesting situation: the kinematical backreaction is 
responsible for two accelerated phases. Our model only concerns the matter-dominated 
universe; however this situation allows to imagine that the primordial inflation and the 
one occuring today might be driven by the same effect studied in a more general model. 

6. Concluding Remarks 

We have built in this paper a model in which the backreaction fluid acts as a domain- 
dependent effective CG, focusing for simplicity on the case where its effective energy 
density is positive. The behaviour of the kinematical backreaction only depends on the 
initial conditions of the domain and may correspond to dark matter or different types 
of dark energy according to the scale and the time evolution. Our model kinematically 
resembles a Friedmannian cosmology with two fluids (dust and CG), but conceptual 
implications differ. First, the origin of dark energy (or, in the extreme case, of both dark 
components) is related to the non-trivial geometrical structure of an inhomogeneous 
space; we do not assume the existence of other fundamental constituents or fields. 
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Second, the geometrical structure that complies with the Chaplygin equation of state 
can in principle be verified by concrete inhomogeneous models and also observationally; 
the model is no longer phenomenological in the sense that no free parameters remain. 
Free parameters in the standard approach (being homogeneous) are here traced back 
to the initial data for the inhomogeneities and are therefore, in this sense, not 
arbitrary. Any fitting of observational data of our model will lead to unambiguous 
initial data that can be constraint by structure formation. Third, one has to be careful 
when cosmological observables are derived. Indeed, angular diameter and luminosity 
distances, for instance, depend on metrical properties. These latter are affected and 
related, in our description, to the averaged scalar curvature, which evolves differently 
compared with its Friedmannian counterpart. It is therefore necessary to reinterprete 
observational data, using e.g. the lines of the analysis performed for the exact scaling 
solutions in [38] . In light of this remark it is premature to exclude [11161171121] the standard 
Chaplygin equation of state as providing a good match with observational data. The 
construction of an effective metric for a cosmology with a Chaplygin backreaction fluid 
and a comparison with observations are the subjects of a future work. 
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Appendix A. Study of the model in the case = 

Setting Vt^ to allows to compare our model to the concordance model (we also suppose 
here that fl^ ^ for any domain). In this situation equation ( 1491) becomes 

n£ + n5, + ng = i nj^i-n,,. (A.i) 

Since and fty^ + Qq are both positive under the Chaplygin fluid constraints (see 
equation (1501)). one also needs 

n£<i, n5, + ng<i. (A.2) 

A.I. Evolution of the Hubble functional 
Equation (|55|) becomes 

1 - n a ) a v 3 + (Q a Qp + 2 tt a tt~, 0^ ) 1/2 j . (A.3) 



Tj2 tt2 



For fi 7 > H%, always decreases, and for fi 7 < it increases in the interval (a™ n , ai) 
and decreases in (ai, +00) with 

<«S*)' = -2g, W « = -2| + 4 ?1 -^L i . (A.4) 
A.2. Evolution of andVLq 

The derivatives of equations (loll) and (j521) show that, for f2 7 > 0, increases and Qq 
decreases. For f2 7 < 0, decreases in (a™ n , 02) and increases in (02, +00), and Qq 
increases in (a™ n ,a 3 ) and decreases in (03, +00), with 
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-Oy [3(1 - tt a ) 2 - 2tt a tt 7 ] + 2Jtt a tt* - (1 - 0*) 2 ] 

«*>' o^fe ' (A ' 5) 

M ° = a^OT ' (A ' 6) 

A. 3. Evolution of the deceleration parameter 
We define the quantities 

(as)6 . = 8 n a n 7 + (i - ap 2 - (i - y/(i - ap 2 + 48 gg; g 

8 fiaf^g 

(ag)6 . = 8 n a n 7 + (i - a,) 2 + (i - n a ) y/(i - a,) 2 + 48 n a n 7 (A 

8 

12 + 1 - 2^6(3^-1) (2^ + 1) 
_^ 12 + 1 + 2y/6(3fy-l)(2fy + l) 

a 2 25 , I • ) 

0*3 := \ - (A.12) 
and the situations 

(a) Va v e « n ,+oo) g© < ; 

(b) Vap G (a™ n , +oo)\{a 4 } qr> < 0, and qv(o-4> = ; 

(c) Vac G (a£> im ,a 5 ) U (a 6 ,+oo) g© < , g©(a 5 ) = g©(a 6 ) = 0, and Va p G 
(a 5 ,a 6 ) > 0; 

(d) Va© G (0, a 6 ) g© > , qv(a 6 ) = 0, and Vap G (a 6 , +oo) g© < . 

Table 1 presents the exact evolution of the deceleration parameter in these 
situations. 



Chaplygin Gas and Inhomogeneous universe 



21 





(0,%) 




0<n (3 <l/3 


(c) 


(d) 





(o,n a3 ) 






[fy,l) 


0.(3=1/3 


(c) 


(b) 


(c) 


(d) 





(o,n ai ) 








(S7 a2 , Q^) 


[%,1) 


l/3<n (3 <l 


(c) 


(b) 


(a) 


(b) 


(c) 


(d) 





(o,n ai ) 




(n Q1 ,i) 


1<% 


(c) 


(b) 


(a) 



Table 1. Evolution of the deceleration parameter for the different values of Q a and 
flp under the Chaplygin fluid constraints with Sl^ = 0. In situations (a) and (b) the 
expansion of the domain is always accelerated; in situation (c) it is first accelerated, 
then decelerated and again accelerated; in situation (d) it is first decelerated, then 
accelerated. 



